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Abstract 

This talk is based on results obtained for masses and wave functions of heavy 
quarkonia in a light-front Hamiltonian formulation of QCD with just one flavor 
of quarks using an ansatz for the mass-gap for gluons. Since the calculated 
spectra compare reasonably well with data, some further steps one can make 
are discussed. 



1 Introduction 

Discussion of heavy quarkonium dynamics in this talk is based on results for 
masses and wave functions obtained in Refs. D 12) in a relativistic (boost- 
invariant) Hamiltonian formulation of QCD. Steps involved in the calculation, 
starting with the Lagrangian for QCD, deriving the corresponding canonical 



light-front (LF) Hamiltonian, carrying out the renormahzation group proce- 
dure for effective particles (RGPEP) to obtain the quark and gluon operators 
at finite momentum scales, A, deriving an effective Hamiltonian for heavy con- 
stituent quarks and gluons, at momentum scales on the order of the quark 
mass, Ao ~ m, using an ansatz for the mass-gap for gluons, to finesse a 
new Hamiltonian that acts only in the effective quark-antiquark Fock sector, 
H^^QQ, and solving numerically the resulting eigenvalue problem for H^^^qq, 
are described in the original literature (a condensed summary 1^1 is available). 

Here, only one example of results for quarkonium masses obtained from 
HxfjQQ is quoted, to illustrate what happens in the simplest version of the LF 
approach to QCD. The key point is that the results do not depend on the ansatz 
^ and fit data reasonably well for the coupling constant expected from RGPEP, 
assuming it has a known value at A = Mz, and for the charm or bottom quark 
masses that have typically considered sizes. Then, the emerging recipe for the 
mass gap ansatz as a tool to facilitate numerical studies of effective quark and 
gluon dynamics is described. The ansatz is designed to be introduced only in 
the final stage of diagonalizing H^^qq. 

LF quantum field theory has a long history ® with lots of modern de- 
velopments ^ that cannot be duly reviewed here. As part of the progress, a 
conceptual outline of nonperturbative QCD in the LF frame was achieved ® 
using the similarity renormahzation group procedure 13 IS! . A confining loga- 
rithmic potential of order a in quark-antiquark sector has been discovered 
and studied in heavy quarkonia and other systems , but not 

using RGPEP ^ . Besides the work that led to Refs. E © , distinct RGPEP 
applications in scattering theory and gluonium are relatively new. But 
there is a lot of work to do in LF QCD before it can be widely accepted as a 
viable alternative to lattice QCD 113). The AdS/CFT method [ISl \M has not 
been connected with RGPEP in QCD yet. 

2 Bottomonium masses as example 

In the crudest version, H\ is calculated using RGPEP to first order in a\g 
in one flavor QCD, and H^^qq is then evaluated to the same order, using 
the ansatz for gluons in the effective quark-antiquark-gluon sector (see next 
section). The resulting eigenvalue problem does not depend on the ansatz ^ 
and takes a form that satisfies requirements of rotational symmetry despite 
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that the LF reference frame distinguishes z-axis. For example, the oigcnstates 
with quantum numbers of T are described by a 2 x 2 matrix wave function 
(p{k ) = b{k )(?, where (s is the polarization three-vector that determines the 
polarization state of the whole quarkonium in motion with arbitrary velocity) 
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The matrix 4) enters into the definition of a relativistic quantum state of a 
quarkonium with definite momentum and mass M, 

\M, P+, P^, s) = J [ij] {2nfP+S'{P -Pi- pj) x\ mj)Xi blAoJ |0) ' 

(2) 

where i and j denote flavor, momentum, and spin quantum numbers (colors 
are combined to 0) of scale-dependent effective quarks that are created by op- 
erators b\^^^ and d\^j from the LF QCD vacuum. These LF operators depend 
on the scale A (in ratio to Aqcd in the RGPEP scheme, and the quark mass). 
The relative momentum three- vector kij is defined using LF kinetic momentum 
variables and the quark mass corresponding to the scale Aq. The eigenvalue 
equation satisfied by the S and D wave functions is written in terms of dimen- 
sionless momentum variables 



p = kij/ki 



(3) 



where ks denotes the strong Bohr momentum, a\^^m,\„/2 (subscript Aq will be 
dropped from now on). Using p = |p|, the radial equation can be written as 
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while the quarkonium mass eigenvalue is given by the eigenvalue x through 



M = 2m 
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Note that the eigenvalue is not energy in any specific frame of reference but the 
mass itself. The functions Wss, VVsd, yVds, V^dd are given in the literature'^. 

There is no quantitative trace of the gluon mass ansatz in this result. 
But there is a qualitatively new element in the form of a harmonic oscillator 
correction to the strong Coulomb potential (with LF Breit-Fermi terms). 

Another qualitatively new element, a result of using RGPEP, is the form 
factor 



/ = exp - 



X2 
-^0 



(8) 



where M denotes an invariant mass of a pair of free quarks. The form fac- 
tor tempers the spin-dependent gluon exchange interaction. In particular, it 
regulates otherwise ultraviolet-divergent three-dimensional delta functions (in 
the position space formally associated with the momentum space of k via the 
Fourier transform), which are present in the functions W due to the relativistic 
spin effects. 

One solves the eigenvalue equations for bb bound states, such as eqdl 
assuming that a is given by the RGPEP evolution from the known value at 
A = Mz down to Aq. If auz = 0.12, the lowest order RGPEP evolution of a\ 
in QCD with only one flavor!^ produces a ~ 0.326 at Aq ^ 3.7 GeV (about 
30% smaller value is generated for 6 or 5 flavors). Less is known about the 
RGPEP evolution and value of the 6-quark mass, mi,. TablT] shows masses of 
bb quarkonia obtained when a and mt, are adjusted to reproduce masses of 
Xi(lP) and Xi(2P) at Aq = 3697.67 MeV. 

If the RGPEP calculation of H\ and subsequent reduction to H^qq were 
exact, there should be no dependence of the spectrum on A. Once a and 
TTib are adjusted to observables at one scale, they evolve in some exact way, 
including the formation of bound states. But in this crudely simplified version 
of LF QCD, the RGPEP procedure is limited to order a and H^^qq is finessed 
using an ansatz for the gluon mass gap /i. Therefore, one cannot change A 
considerably using equations limited to order a and there is only a hope that 
in some small range of values of Aq the equations have a chance to work once the 
coupling constant and quark mass are given their right relativistic values^ . 

There is a characteristic pattern visible in the fourth column in tab[TJ the 
greater the difference between a mass eigenvalue and the masses of quarkonia in 
the middle of the table, used to choose a and rrih, the greater the discrepancy 
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Table 1: Example of calculated masses (MeV) for bb states. The corresponding 
coupling constant and quark mass are a = 0.32595 and mi, = 4856.92 MeV. 



meson 


theory 


experiment 


difference 


T10860 


10725 


10865 ±8 


-140 


T10580 


10464 


10579.4 ±1.2 


-116 


T3S 


10382 


10355.2 ±0.5 


27 


X22P 


10276 


10268.65 ±0.22±0.50 


7 


Xi2P 


10256 


10255.46 ±0.22±0.50 





Xo2P 


10226 


10232.5 ±0.4 ±0.5 


-6 


T2S 


10012 


10023.26 ±0.31 


-11 


X2IP 


9912 


9912.21 ±0.26±0.31 


-1 


XilP 


9893 


9892.78 ±0.26±0.31 





XolP 


9865 


9859.44 ±0.42±0.31 


5 


TIS 


9551 


9460.30 ±0.26 


91 


77f,lS 


9510 


9300 ±20 ±20 


210 



between the crudely approximated theory and experiment. This should be 
expected. The most strongly bound states are sensitive to deviations of the 
effective potential from the Coulomb shape. For example, interactions order 5^ 
(or a^), introduce ^- functions that are absent here because of the limitation to 
terms order in the RGPEP and two-quark reduction. Analogous 4th order 
6s and other singular corrections are known in QED. Here such terms should 
have much larger effect because the coupling constant is about 30 times larger 
then in QED (they have to be treated nonperturbatively). The least strongly 
bound states, those with largest masses, should not be described well without 
proper inclusion of gluons. The mass ansatz should fail to render interactions 
that are associated with gluons producing a linear potential at distances much 
larger than the strong Bohr radius. 

3 Mass ansatz as a computational tool in QCD 

The mass ansatz for virtual effective particles in Fock components that contain 
more such particles than just two in mesons, or three in baryons, deserves a 
comment for several reasons. One reason is that the ansatz may help improve 
the calculations for heavy quarkonia. Another reason is that it may lead to a 
possibility of calculating properties of baryons built from heavy quarks. The 
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third reason is that it may help in crossing the barrier that separates ah smaU- 
couphng expansions in QCD from entering the region of quark masses much 
smaller than Aqcd- 

The first two reasons concern the difficulty that precise numerical solu- 
tions of eigenvalue problems with coupled two-, three-, and more-particle sec- 
tors are hard to obtain. In this respect, the mass ansatz appears a candidate 
to mimic what happens in an infinite tower of the Fock components built from 
effective particles. Since the RGPEP form factors limit momentum transfers 
by A, the spread of probability to sectors with many effective particles corre- 
sponding to the scale A is tamed. But these sectors do infiuence the dynamics 
of the dominant sectors and some ansatz appears inevitable. The question is 
how to make it self-consistently. The basic idea is to drop all sectors above 
the highest included (in the sense of number of the effective particles), put in 
instead a mass ansatz in the highest Fock component, and see what happens 
in the dynamics of lower components. The next step is to increase the maxi- 
mal number of effective particles by one and see if the same type of ansatz is 
producing the same answers. A small coupling expansion may constrain the 
options sufficiently for finding good candidates for suitable mass terms in the 
highest components. 

The third reason that concerns light quarks is most speculative. It in- 
volves chiral symmetry, or rather the mechanism of its breaking. In LF Hamil- 
tonian of QCD with small A, there may exist finite terms that violate chiral 
symmetry and do not vanish in the limit of quark mass approaching zero ^ . 
At the same time, the LF vacuum state remains simple due to cutoffs imposed 
on the particle momentum components along the front. The question is how to 
find those terms in practice. LF power counting limits the structure of allowed 
terms but so far insufficiently for anybody to tackle the issue, even though the 
stakes are high. 

From this point of view, the following observation is of interest. Consider 
a colorless state built from two effective gluons. They attract each other. 
Consider then that one of these gluons turns into a pair of quarks. These 
quarks are in an octet state and instead of attracting they repel each other. 
In perturbation theory, if the number of quark flavors is not too high, this is 
not dangerous and gluonic interactions sustain asymptotic freedom, generating 
infrared slavery. However, beyond perturbation theory, a pointlike creation of 
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a pair of quarks that repel eaeh other by violent potentials may lead to an 
explosive behavior. Such behavior is entirely absent in QED because electrons 
attract anti-electrons and this effect slows down the growth of a pair, instead 
of accelerating it. 

To be more specific, consider eqlH in which the s-wave potential Wss 
contains a term © that would be a (5-function if the form factor / did not 
smooth it out. This term is attractive in T. It is repulsive instead in the color 
octet states. When A is comparable to the heavy quark mass, which means 
that a is small, the smallness of the coupling constant and form factor width 
produce an interaction that cannot compete significantly with the size of the 
quark mass. This is visible in tab[TJ However, if the quarks are light, it is 
entirely unclear what will happen. 

The situation is different than in the case of analogy between a gluonium 
and a helium atom with one doubly charged electron discussed in Ref. ^^21) _ 
Here, two particles with the same charge are suddenly put on top of each other 
and large potential energy is created. 

In the case of light quarks, the large terms are smoothed out by the 
RGPEP form factors / but their strength may be comparable with Aqcd- 
The central point is that in order to find out what happens due the explosive 
nature of color dynamics of effective particles beyond perturbation theory, one 
has to separate some sectors from a presumably decreasing but in principle 
infinite chain of them. This is what the mass ansatz facilitates. Thus, it opens 
a way to investigate interactions in the effective Hamiltonians that may be 
calculated perturbatively in RGPEP and then diagonalized nonperturbatively 
using computers. This way one can find out if the effective interactions may 
in principle be responsible for emergence of the constituent quark masses for 
quarks u, d, and s. 

4 Conclusion 

Tab[T] shows that the chance that LF Hamiltonian approach to QCD may 
apply in phenomenology of heavy quarkonia is not hopelessly small. Since 
the coupling constant one needs is order 1/3, and one may need even a smaller 
coupling constant when 4th order RGPEP is used, the approach stands ready 
for a more extended scrutiny. The reason it deserves to be checked is that 
it appears now to indicate a possibility that a single formulation of the entire 
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theory with quark masses much greater than Kqcd is conceivable with no need 
to combine different formulations for including information concerning different 
scales. 

The harmonic potential finessed in the quark-antiquark sector using an 

ansatz for the gluon mass gap, leads to the eigenvalues that are propor- 
tional to the angular momentum of relative motion of quarks, like in the Rcgge 
trajectories. It is found that the oscillator frequencies are on the order of one 
inverse fermi, and the oscillator potential grows as the relative distance squared 
in fermis with a coefficient given by the quark mass. So, presumably, for states 
with masses greater than order 1 GeV above the ground states, the probability 
of emission of effective gluons increases and then formation of strings of gluons 
is favored if the gluons also have some oscillator force acting among them. For 
quantum gluons to form a string, each pair of the neighboring gluons must be 
held together stronger than by a linear potential. If it is capable of generating 
such effects, the LF Hamiltonian approach could thus lead to a quantum theory 
of the gluon string in QCD without ever introducing a nontrivial vacuum. 

But the concept of mass ansatz in Fock sectors with one more effective 
particle than the maximal number treated nonperturbatively, is probably most 
interesting as a tool for finding out what happens when one attempts to solve 
eigenvalue equations for Hamiltonians that are evaluated using perturbative 
RGPEP in the case of canonical QCD with quark masses smaller than Aqcd- 
The idea discussed here is that effective particles in non-singlet color configura- 
tions may experience explosive potentials in the form of smoothed 5-functions 
that may cause effects order Aqcd per constituent. The ansatz for quark and 
gluon masses inserted in the highest sectors thus opens a possibility to generate 
concrete forms of such terms and to study them nonperturbatively. Nothing is 
known yet about what may come out from such studies. 
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